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A BS TRACT 

Recent discoveries in material science and 
fluidics have been used to create a variety of 
novel effector devices that offer great potential to 
enable new approaches to aerospace vehicle 
flight control. Examples include small inflatable 
blisters, shape-memory alloy diaphragms, and 
piezoelectric patches that may be used to pro- 
duce distortions or bumps on the surface of an 
airfoil to generate control moments. Small jets 
have also been used to produce a virtual shape- 
change through fluidic means by creating a recir- 
culation bubble on the surface of an airfoil. An 
advanced aerospace vehicle might use distributed 
arrays of hundreds of such devices to generate 
moments for stabilization and maneuver control, 
either augmenting or replacing conventional 
ailerons, flaps or rudders. This research demon- 
strates the design and use of shape-change device 
arrays for a tailless aircraft in a low-rate maneu- 
vering application. A methodology for assessing 
the control authority of the device arrays is 
described, and a suite of arrays is used in a 
dynamic simulation to illustrate allocation and 
deployment methodologies. Although the 
authority of the preliminary shape-change array 
designs studied in this paper appeared quite low, 
the simulation results indicate that the effector 
suite possessed sufficient authority to stabilize 
and maneuver the vehicle in mild turbulence. 


INTRODUCTION 

A number of concepts for novel control 
effector devices have recently been developed 
which offer great potential to enable entirely new 


approaches to aerospace vehicle flight control. 
Examples of such devices include small 
inflatable blisters, shape-memory alloy 
diaphragms and piezoelectric patches that may 
be used to produce distortions or bumps on the 
surface of an airfoil, as well as small 
piezoelectric pumps that can produce oscillatory 
surface jets. Potential applications of these 
devices include forebody vortex asymmetry 
management (Roos), control and manipulation 
of shear flows (Smith, Glezer), fluidic thrust 
vectoring (Pack), and low angle-of-attack lift 
augmentation (Nae). 1 ' 4 Pulsed blowing has also 
been used to actively vary the reattachment 
location of a separated flow over specifically 
designed portions of an airfoil (Seifert, Pack). 5 6 
The Aircraft Morphing Program under way at the 
NASA Langley Research Center focuses on the 
development and application of these novel 
effector devices, covering a spectrum research 
topics which include materials development, 
actuator design and testing, CFD modeling and 
optimization methods, dynamic simulation and 
control using such devices (McGowan). 7 In this 
study we will examine the potential to use 
distributed arrays of shape-change devices to 
provide a basis for aircraft flight control. An 
advanced aerospace vehicle might use hundreds 
of small shape-change devices to generate forces 
and moments for stabilization and maneuver 
control, without the need for conventional, 
hydraulically actuated ailerons, flaps or rudders 
as investigated by Scott and Montgomery, and 
by Allan. 8,9 Potential reasons to consider such 
alternatives include reduced observability, 
redundancy, and possible weight reduction of the 
overall effector complement. 
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Given this goal, how might the best 
locations for placement of such devices on the 
aircraft surface be determined, and how might 
the devices be used in a control system to 
actually maneuver and fly the vehicle? This 
research attempts to address such questions by 
applying arrays of a generic shape-change device 
to a representative aircraft configuration in a 
dynamic simulation. This work builds upon an 
initial investigation that was conducted by Scott 
and Montgomery 8 , extending the research to 
include an effector placement study followed by 
a revised effector suite, a new allocation and 
deployment method, and a new control design. 
The control system deploys distributed arrays of 
the shape-change devices in a “quantized” 
fashion. That is to say that each small bump in 
an array is either completely on or off, and more 
of them are turned on to produce larger forces as 
needed. Using these arrays, the control system is 
able to stabilize and maneuver the vehicle 
without conventional hinged surfaces such as 
ailerons or a rudder. The predicted authority of 
the designs studied in this paper is rather low 
when compared to conventional effectors, so the 
target for utilization of this technology is 
currently low-rate maneuvers (roll rates of 10 
degrees per second or less). 

SUPJE CT CO NF IGURATION 

The example aircraft design used in this 
paper is the Innovative Control Effector (ICE) 
aircraft from Lockheed Martin Tactical Aircraft 
Systems, (Dorsett). 10 The configuration is a 
tailless, delta-wing fighter, which was developed 
for the purpose of investigating novel control 
effector concepts and arrangements. NASA is 
using the ICE design, shown in Figure 1, as an 
example configuration under a cooperative 
agreement with Lockheed Martin. 

Wing Characteristics 

Area .... 75.12 m 2 (808.6 ft. 2 ) 

Span .... 11.43 m (37.5 ft.) 

Aspect Ration ... 1.74 

Leading Edge Sweep .. 1.134 rad 


Figure 1. Lockheed- Mart in's ICE configuration . 

A simulation database of stability and 
control derivatives was obtained for this 
configuration from Lockheed Martin. The 
potential of shape change device arrays to serve 
as alternative control effectors for the ICE 



aircraft is examined using the simulation model 
with linearized aerodynamics. A low-speed, 
potential flow code (the Panel Method from 
Ames Research Center, PMARC) is used to 
predict the forces and moments that would be 
produced by the shape-change effector arrays 
which were designed for this configuration. 11 

The analysis was performed for a flight 
condition of Mach 0.6 at an altitude of 15000 ft. 
The equilibrium trim angle-of-attack at this flight 
condition is approximately 4.4 degrees. Control 
moments required to trim at this flight condition 
are assumed to be generated by conventional 
effectors. The shape-change effector arrays are 
used to stabilize and maneuver the vehicle about 
this equilibrium flight condition. In this paper, 
the term "device” will be used to refer to an 
individual shape-change effector unit or single 
bump, the term "array" will refer to a particular 
grouping of adjacent devices, and the term 
"suite" will refer the entire collection of device 
arrays with which the vehicle is outfitted. 

CONTROL EFFECTOR DEVICE MODEL 

The generic shape-change device model 
used in this investigation consists of a small 
bump which deflects normal to the surface of the 
vehicle. The resulting shape-change was 
represented on a panel model of the ICE 
configuration by displacing a given gridpoint a 
prescribed distance, X N , along a vector which is 
normal to the geometry surface, as shown in 
Figure 2. The perturbed geometry was then 
analyzed using the PMARC code, and the 
resulting forces and moments were differenced 
from the values for the nominal geometry to 
provide an estimate of the forces and moments 
that would be produced by the surface distortion. 
Several assumptions are implicit in this modeling 
approach, including the assumption that the 
surface distortion does not produce separated 
flow and the assumption that the flow is entirely 
subsonic. These restrictions were imposed by 
the potential flow code that was used to evaluate 
the shape-change devices. 



Figure 2. Shape-change device modeled as 
deflection of gridpoint along normal vector. 
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The use of the modeling approach illustrated 
in Figure 2 appears reasonable for shape- 
memory alloy diaphragms, piezoelectric domes 
and inflatable blisters which produce real surface 
distortions. The use of fluidic devices, such as 
pulsed jets and zero-net-mass-flow oscillatory 
jets ("synthetic jet"), to produce a virtual shape- 
change is far less understood and is a subject of 
current research. Krai and Donovan have 
performed a detailed study regarding the 
simulation of an isolated jet using Reynolds- 
averaged Navier-Stokes (RANS) equations. 12,13 
Their result, shown in Figure 3, demonstrates the 
effect of a steady-jet on a NACA 0012 airfoil at 
r n = 8.5 x 10 6 and a =.0698 rad (4.0°), The 
figure shows streamlines for the U/LL = 1.0 
case, illustrating that the effect of the steady jet 
is to produce a virtual distortion of the surface 
resulting from the generation of a small 
recirculation bubble, as well as an effluent added 
to the primary flow. A proper potential model 
would include a surface distortion as well as an 
effluent. In this paper we assume there is zero- 
mass-transfer and therefore no added effluent. 



Figure 3. Streamlines a= .0698 rad (4.0 °), with 
steady-jet control at Uj/U„ = 1.0. 

The modeling approach described above 
includes no description of the actuator dynamics. 
A recent investigation by Pack found the 
response time of an exhaust flow to a thrust- 
vectoring pulsed jet to be very fast -- on the 
order of 22 ms, including both the actuation of 
the jet device itself and the fluidic response to 
the presence of the jet. 3 Other means of 
generating the shape-change, such as inflatable 
blisters, could presumably be slower since the 
response time would clearly be dependent upon 
the particular mechanization employed. The goal 
of this investigation was to examine the issues of 
shape-change effector placement, deployment, 
and use in control design, rather than to develop 
a high-fidelity model of a shape-change effector, 
and so an instantaneous response of the device 
was deemed sufficient for this first-order model. 


DESIGN OF EFFECTOR ARRAYS 

A major issue to consider in the design of 
shape-change effector arrays is where to locate 
the devices on the aircraft surface for maximum 
benefit. A sensitivity analysis was performed 
using a geometry model of the ICE configuration 
to aid in determining where the device arrays 
should be positioned on the vehicle. The analysis 
involved differentiating the PMARC code using 
the "ADIFOR" tool (Automatic Differentiation 
of Fortran, Carle & Fagan). 14 The differentiated 
code was then applied to the ICE configuration 
model to produce the partial derivative of the 
forces and moments acting on the aircraft with 
respect to a displacement along the surface 
normal at every grid point on the aircraft 
geometry model, as described by Park & 
Green. 15 The resulting sensitivity database of 
pitch, roll, and yaw moment derivatives with 
respect to a normal displacement at every 
gridpoint on the geometry (taken individually) is 
shown as a set of shaded planform maps in 
Figure 4. The shaded regions in the figure 
indicate greater moment sensitivity and therefore 
suggest potential locations at which to place the 
shape-change devices. The sensitivity plots are 
described in much greater detail by Park & 
Green. 15 




Figure 4. Moment sensitivity databases for 
gridpoint displacements produced by applying 
ADIFOR to PMARC using the ICE model 
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The databases shown in Figure 4 were 
incorporated into an interactive design tool that 
was developed to allow an engineer to quickly 
build up and analyze distributed arrays of small 
shape-change devices on the surface of the ICE 
aircraft. The tool was produced using Matlab™ , 
which provides basic constructs for generating 
the interactive graphical user interface and data- 
manipulation functions. It allows the designer to 
use the sensitivity data in a subjective fashion to 
guide the placement, size, and shape of the array 
so that it can produce the forces and moments 
required to maneuver the vehicle. 

The graphical user interface for the effector 
array design tool is shown in Figure 5. Using 
this tool, the designer designates the geometry 
gridpoints at which to place each element of the 
shape-change array. The designer also defines a 
unit deflection magnitude for each element of the 
array. Once a grouping of shape-change devices 
has been defined, the designer can obtain a 
preliminary prediction of its effectiveness and 
generate a perturbed geometry grid that includes 
the deployed effector array. This geometry file 
can then be used with aerodynamic analysis 
programs, in this case the PMARC code, to 


further assess the effectiveness of the device. 
Although the design tool provides some 
subjective guidance regarding the location and 
shape of a given device array, it does not provide 
insight regarding the appropriate height profile 
(deflection magnitudes) for the various elements 
of the array. Such insight might be provided by 
a matrix of second derivatives of the control 
moment with regard to normal displacement of 
the geometry gridpoints (the Hessian). The 
ADIFOR automatic Fortran differentiation code 
has recently been endowed with a capability to 
generate the Hessian of the dependent 
parameters in a Fortran program by Carle and 
Fagan of Rice University, but this new capability 
has not yet been applied to the ICE shape-change 
sensitivity analysis. 

The design tool was used to generate a series 
of 31 candidate effector array configurations. A 
plot of the pitch, roll, and yaw effectiveness for 
each of the arrays is shown in Figure 6. 
Candidate arrays were designed for various 
regions on the configuration, including upper 
and lower surface leading and trailing edges, 
mid-chord locations, and several wingtip 
positions. 



Figure 5. Shape -change effector array design tool graphical user interface . 
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Figure 6. Control authority plots for 31 candidate shape-change effector array designs . 
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The candidate designs were evaluated on the 
basis of their authority and their ability to 
generate a relatively decoupled control moment. 
The total number of devices included in a given 
array was not considered as a factor in the 
evaluation. A particular combination of four 
device arrays was selected based on a subjective 
evaluation of suitability for providing 3-axis 
moment control. The effectiveness derivatives 
corresponding to the selected array designs are 
circled in Figure 6. An ongoing investigation 
being performed in the Multidisciplinary 
Optimization Branch at the NASA Langley 
Research Center is examining the application of 
a more rigorous optimization method for the 
selection of an effector array suite from the 
candidate designs. 16 The suite of four effector 
arrays that was selected for use in the remainder 
of this study is described below. 

SELECTED EFFECTOR ARRAYS AND 
DEPLOYMENT SCHEME 

The four distributed shape-change device arrays 
that were selected to compose the effector suite 
for the ICE configuration are shown in Figure 7. 
The effector suite includes four arrays on each 
wing, three on the upper surface and one on the 
lower surface. The entire suite of shape-change 
effector arrays includes a total of 156 individual 
devices, 78 per wing. The upper-surface leading- 
edge array (ULE) consists of 10 devices - one 
device at each of 10 span stations (chord lines) 


distributed along the leading edge as shown in 
Figure 7. The lower-surface trailing-edge array 
(LTE) and upper-surface trailing-edge array 
(UTE) both consist of 22 devices - two adjacent 
devices at each of 1 1 span stations. The upper- 
surface wingtip array (UTip) consists of 24 
devices - six adjacent devices at each of 4 span 
stations. The lower-surface trailing-edge, upper- 
surface trailing-edge, and upper-surface leading- 
edge arrays all have a maximum displacement 
height of 2.4 inches. The upper-surface wingtip 
array has a maximum displacement height of 1.2 
inches. 


Upper Surface 
Trailing-Edge 
(UTE) Array 


Upper Surface 
Wingtip (UTip) 
Array ^ 



Upper Surface 
Leading-Edge 
(ULE) Array 


Lower Surface 
Trailing- Edge 
(LTE) Array 


devices activate 
inboard to 
outboard as 
more controf is 
required 



Figure 7. Selected shape-change effector array 
designs applied to ICE configuration. 
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The effector suite was applied to the ICE 
vehicle in a simulation and used in a stability 
augmentation and control system design. The 
control system employs a “quantized spanwise” 
deployment scheme for activation of the device 
arrays, as noted in Figure 7. That is to say that 
each device in an array is either completely on or 
off (rather than proportionally deployed), and 
more of them are turned on to produce larger 
forces as needed, from the most inboard to the 
most outboard span station of each array. All 
devices located along a particular span station in 
a given array operate together - turning on or off 
at the same time. This quantized deployment 
scheme addresses a realistic aspect of the control 
allocation problem, since a number of existing 
shape-change device concepts involve this sort 
of "bang-bang" operation. 

Each of the four arrays were analyzed using 
PMARC in each of their quantized deployment 
states to generate a database of control 
effectiveness for use in a closed-loop dynamic 
simulation of the vehicle model. Figure 8 shows 
plots of the control effectiveness for the 
progressive spanwise deployment of each of the 
device arrays. The use of the PMARC 
predictions for control effectiveness in the 
dynamic simulation constitutes a quasi-static 
assumption since unsteady aerodynamic effects 
due to the transient deployment of the devices 
are neglected. 

p yna mic simulatio n 

A six degree-of-freedom dynamics 
simulation was used to investigate the 
unaugmented and augmented aircraft dynamics. 
The body axis coordinate system shown in 
Figure 9 is used for the formulation of the 
equations of motion. 

The longitudinal configuration aero- 
dynamics are represented in the simulation as a 
linear expansion of the longitudinal aerodynamic 
force and moment coefficients in angle of attack, 
a, and pitch rate, q, plus the force and moment 
contributions from the shape-change effector 
deflections, 5, as in the method of Bryan. 17 Thus, 
the longitudinal coefficients, Cx, Cz, and Cm, 
are expressed as: 

Cx = Cx 0 + Cx a oc + Cxqq(c/2 V) + Cx 6 5 
Cz = Cz 0 + Cz a a + Cz q q(c/2V) + Cz 5 8 
Cm = Crrio + Cm a a + Cm q q(c/2V) + Cm 6 8 






Number of Span Stations Activated 

Figure 8. Control effectiveness buildup for each 
of 4 shape-change effector arrays included in 
ICE effector suite (right wing). 
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X, p 

Cx, Cl 


Figure 9. Body axis coordinate system 
definitions used in dynamic simulation . 

Dependence on both angle of attack, a, and 
angle of sideslip, P, is retained for the lateral- 
directional force and moment calculation, so the 
coefficients, Cy, Cl, and Cn, are expressed as: 

Cy = Cy p P + Cy pP (b/2 V) + Cy r r(b/2 V) + Cy 5 8 

Cl = CluP + Cl pP (b/2V) + Cl r r(b/2V) + Cl§5 

Cn = Cn p p + Cn pP (b/2V) + Cn r r(b/2 V) + Cn 6 5 

Where, 

Cyp = Cy po + Cy^a 
Clp = Cl^ 0 + ClpdOt 

Cn p = Cnp 0 + Cnj^ct 

The parameters relating to flight condition, 
including total air speed, V, and dynamic 
pressure, qbar, are calculated as: 

V = (u 1 2 + v 2 + w 2 ) 1/2 
qbar = p*V 2 /2 
a = arctan(w/u) 

P = arcsin(v/V) 


The geometry constants, c, b, and S are the 
reference chord length, wing span, and wing 
area, respectively. Also (u,v,w) and (p,q,r) are 
the translational and rotational rates about the 
body-fixed, reference (x,y,z) axes of the aircraft, 
respectively, and (Cx,Cy,Cz) and (Cl,Cm,Cn) are 
the aerodynamics force and moment coefficients 
about the (x,y,z) axes, as indicated in Figure 9. 
The equations of motion for the rigid body 
dynamics of aircraft are well developed and 
available in texts (Etkin). 18 They require a 
description of the attitude of the vehicle, which 
has traditionally been via Euler angles. These are 
used herein to represent the orientation of the 
aircraft with the sequence corresponding to \| f 
about the z the axis, 0 about the y axis, and 
<j) about the x axis. The kinematic differential 
equations implemented in the simulation are: 


d0/dt = q*cos(<|>) - r*sin(<J>) 

dy/dt = (q»sin(<()) + r«cos(<()))/cos(0) 

d<|>/dt = d\|//dt*sin(0) + p 

dh/dt= u , sin(0)-w # cos(0)cos((t>)-v , cos(9)sin((|)) 

where h is the altitude. For the longitudinal 
variables, the dynamic equations of motion used 
are: 

du/dt = v»r-w*q+g[(T +Cx*qbar*S)/W- sin(0)] 
dw/dt = u # q-v*p+g[cos(0)cos(<|))+Cz'qbar B S/W] 
dq/dt = [Cm*qbar*S*c+(I z -I x )np+(r 2 -p 2 )I xz -r»H T ]/I y 

where g is gravity. Due to inertia terms it is 
convenient to define the total rolling moment, L, 
and yawing moment, N, for the lateral equations 
of motion as: 

L = CFqbar«S*b + (I y - I 2 )*q*r + p«q B I xz 
N = Cn*qbar«S«b + (I x - T y )*q*p - q*r»I xz + q*H T 

where T is the engine thrust, and H T is the 
moment of momentum about the x axis of the 
engine rotors which are assumed to be aligned 
with the body x axis. Assuming symmetry with 
respect to the x-z plane, I x , I y , I 2 are the moments 
of inertia about the x,y,z axes and I xz is the 
relevant product of inertia about the y axis. 
Thus, the lateral equations of motion are: 

dv/dt = w*p-u*r +g(Cy*qbanS/W+cos(0)sin(<!))) 

dp/dt = (I z -L + I X z‘N)/(I x *I z - I X2 2 ) 

dr/dt = (I xz -L + I x -N)/(I x -I z -I xz 2 ) 

The geometry constants used in the 
simulation are c = 28.75 ft., b = 37.5 ft., and S = 
808.6 ft. 2 , and the mass/inertia parameters are W 
= 32750 Ibf, I x = 35479 slug-ft. 2 , I y = 78451 slug- 
ft. 2 , I z = 1 10627 slug-ft. 2 , and I xz = -525 slug-ft 2 . 
The aerodynamic coefficients for the ICE 
configuration at the selected design flight 
condition of Mach 0.6 and 15000 ft are listed in 
Tables 1 and 2. 


Table 1. Longitudinal Aerodynamic Coefficients 
for the ICE simulation. 




Subscript 


0 

a 

q 

Cx 

0.0166 

-0.1973 

0 . 

Cz 

0.0395 

-2.2475 

0 . 

Cm 

0.0036 

-0.0467 

-.39516 
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Table 2. Lateral Aerodynamic Coefficients for 
the ICE simulation. 




Subscript 



(Jo 

Pa 

p 

r 

Cy 

-0.0534 

0.2331 

0. 

0. 

Cl 

0.0109 

-0.7846 

-.016 

.021368 

Cn 

-0.0099 

-0.1215 

-.021789 

-.01 


CO NT ROL DESIG N U SI NG SH APE- 
CHA NG E EF F ECTOR A RRAYS 

Eigenvalues for the uncontrolled flight 
dynamic modes of the ICE aircraft model were 
obtained from a linearization of the simulation 
about the equilibrium trim condition. Table 3 
presents the mode designations, frequencies and 
damping values for the open-loop poles that are 
plotted as circles in Figure 10. 


Table 3. Open-loop poles , frequencies and 
damping ratios for the linearized ICE simulation 
at Mach 0.6 and 15000 ft. 


Eigenvalue 

Damping 

Ratio 

Freq. 

(rad/s) 

Mode 

-8.20e-01± 

3.74e-01 

2.20e+00 

Short- 

2.04e+00i 



Period 

-1.46e-03± 

2.20e-02 

6.64e-02 

Phugoid 

6.64e-02i 




-1.28e+00 

1.00e+00 

1.28e+00 

Roll 

5.72e-01± 

-5.96e-01 

9.60e-01 

Dutch- 

7.71e-01i 



Roll 

-2.27e-02 

1.00e+00 

2.27e-02 

Spiral 



Figure 10. Open-loop and closed-loop pole 
locations for the ICE simulation. 


The longitudinal phugoid and short period 
modes are stable, but the lateral Dutch-roll mode 
is not. The lateral instability is not unexpected 
due to the tailless configuration of the ICE 
aircraft. 

A feedback control law was designed that 
uses the suite of shape-change device arrays to 
stabilize the vehicle's lateral dynamics while 
providing a bank-angle command tracking 
capability. To accomplish the design, it was 
necessary to develop a feedback architecture, 
associated filters and control gains, and to devise 
a control allocation method and deployment 
scheme for the suite of shape-change effector 
arrays. The controller stabilizes and maneuvers 
the vehicle using only the shape-change arrays, 
although the vehicle’s steady-state trim is 
provided by conventional effectors. 

A linear-systems pole-placement approach 
to the control law synthesis was selected for 
simplicity. Since the open-loop longitudinal 
dynamics are stable, the control design focuses 
on the augmentation of the lateral-directional 
dynamics. Therefore, only lateral dynamic states 
were used in the pole-placement controller. An 
overview of the dynamic simulation architecture 
showing the relation of the command generator, 
control law and allocator, is shown in Figure 1 1 . 



Figure 11. Block-diagram overview of 
simulation architecture. 

The linearized lateral dynamics of the open- 
loop aircraft model are represented in state-space 
form as: 


x = [A] x + [B] u 

where the transpose of the state vector, x, is 
given by: 

X T = ( P P r <t> } 

and the transpose of the control deflection 
vector, u, is given by: 

U T “ ( 5uthr SltER $ULER ^UTipR • • ■ 

^UTEL ^lTEL ^ULEL 5 UTipL } 
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where each subscript corresponds to a particular 
effector array as defined in Figure 7, and the 
suffix R or L refers to the right or left wing, 
respectively. The matrices, A and B, contain the 
aircraft stability coefficients and control 
coefficients, respectively, from the linear 
simulation model. To generate the closed-loop 
system we define a commanded input vector 
which drives the shape-change device arrays 
using the following control scheme: 

u cmd = [P] m cmd 

The vector, rn cmd , represents the vector of desired 
roll, pitch, and yaw moments: 

m C md= (Cl Cm Cn) cmd T 

where Cm cmd = 0 and {Cl Cn} cmd is calculated 
by the following lateral-directional state- 
feedback control law: 

{Cl Cn) cmd T = [K] x 

The feedback gain matrix, K, is calculated using 
the pole-placement algorithm for linear systems 
as described by Franklin and Powell. 19 The 
matrix, P, represents the control allocation 
function, which calculates a combined 
deployment of the eight effector arrays that will 
generate an approximation to the desired 
moment vector. In this design, the pseudo- 
inverse allocation method is used, which is 
perhaps the most familiar and fundamental 
means of control allocation for coupled moment- 
generating effectors. 20 The pseudo-inverse 
matrix, [P], of the control effectiveness matrix is 
calculated as: 


P = B t [BB t ] 1 

But since the shape-change effectors are 
one-sided, negative deployments of the devices 
are undefined. Therefore, negative deflections 
are reflected to the corresponding effector array 
on the opposite wing as positive deflections in 
the effector allocation algorithm. An error in 
pitch moment is generated by the reflection of 
negative asymmetric control deflections to the 
opposite wing. This effect is compensated by 
calculating an appropriate symmetric deflection 
of the upper or lower trailing edge arrays, 
depending upon the sign of the pitch moment 
error. 

For the closed-loop system, the following 
desired lateral/directional pole locations were 
prescribed as targets for the pole-placement 
controller: 


Despoles = 


-2.2500 

-7.2500 

-1.7678e + 00 + 1.7678e + 00i 


-1.7678e + 00 - 1.7678e + 00i 


These pole locations, shown as pluses on the 
plot in Figure 10, were selected to stabilize the 
vehicle's lateral dynamics and provide improved 
damping and settling time for gust response. For 
turn coordination, a washout filter is applied to 
the yaw-rate feedback with a transfer function 
that is given by: 


r in 1-5S+1 

A final provision for the use of the effector 
arrays is the presence of limiters within the 
control law, with limits that are dictated by the 
low authority of the control effector arrays. The 
derivative, Cnp, generates an adverse yawing 
moment in proportion to roll rate, tending to 
produce an undesired sideslip transient during 
rolling maneuvers. Additionally, as the vehicle 
rolls about the body axis, angle of attack is 
converted into sideslip. Feedback paths in the 
control law are provided to counteract these 
effects so as to generate coordinated turns. But 
due to the limited yaw authority of the effector 
suite, the ability of the control system to reject 
these effects is minimal, and maintaining turn 
coordination while rolling becomes the limiting 
factor in our achievable roll-rate command. 
Therefore, a roll-rate limiter was placed in the 
bank-angle command path. Iterative simulation 
runs which examined various values for the roll- 
rate limit have revealed that a maximum limit of 
10 degrees per second will prevent loss of 
control during rolling maneuvers, so this value 
was used as the maximum roll rate that could be 
commanded. 

SIMULATION RESULTS 

The effectiveness and utility of the shape- 
change effector array suite was evaluated by 
subjecting the dynamic simulation to various 
atmospheric disturbance levels and by assessing 
the maneuver capability that it provided. Using 
the shape-change arrays, the control system was 
able to stabilize and maneuver the vehicle 
without conventional control surfaces. The 
predicted authority of the shape-change arrays is 
rather low when compared to a rudder or aileron, 
so the control system generates relatively low- 
rate maneuvers. Figure 12 presents two time 
histories produced by the dynamic simulation in 
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response to a i 20-degree bank angle doublet 
command with and without mild Dryden-spectra 
turbulence (a = 3 ft/s). 




Figure 12. Simulation time histories in response 
to a 20-degree bank angle doublet command 
with and without mild turbulence . 


The maximum roll-rate achieved is 10 deg/s 
due to the presence of the limiters in the 
command and feedback paths. Sideslip perturba- 
tions occur during the execution of the roll 
doublet. Larger sideslip transients result from 
the turbulence, but the bank-angle command 
time history is minimally impacted. In spite of 
the relatively low authority of the shape-change 
effector array suite, the control law is able to 
stabilize and maneuver the vehicle within the 
prescribed roll-rate limits. 

Figure 13 shows a comparison of control 
effector array deployment with and without the 
Dryden turbulence. The vertical axis on the plots 
represents the number of devices that were 
activated in each of the eight arrays at any given 


time during the maneuver. The figure is not 
intended to convey specific details regarding any 
one particular array, but rather provides a general 
indication of the level of control activity and 
saturation that occurred during the run. 
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Figure 13. Shape-change effector array 
deployment time histories in response to a 20- 
degree bank angle doublet command with and 
without mild turbulence. 

In both cases, various arrays in the control 
effector suite saturate, that is to say that all 
elements of a given array are turned on, for much 
of the time history. A higher degree of control 
activity and transient saturation results when 
mild turbulence is present, but the control system 
is still able to stabilize the lateral dynamics and 
perform the roll doublet. However, when the 
Dryden turbulence level was increased to 
moderate (<7 = 6 ft/s ), the simulation diverged. 
The moderate turbulence field introduced 
perturbations from the trim condition which the 
effector suite did not have sufficient authority to 
reject, as deployed by the pole-placement control 
law design. 

Time histories produced by the dynamic 
simulation in response to increasing crosswind 
gusts are shown in Figure 14. The vehicle is 
able to recover from crosswind gusts of up to 28 
ft/s, but beyond that the effector suite lacks 
sufficient yaw authority to recover from the 
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bank angle, deg sideslip, deg 


disturbance. This magnitude of crosswind gust 
equates to a sideslip angle of approximately 2.5 
degrees at the simulation flight condition of 
Mach 0.6 and 15 000 ft. The figure shows the 
response to increasing crosswind gust 
magnitudes in 5 ft/s increments, with the 30 ft/s 
crosswind gust resulting in a loss of control. 




Figure 14. Simulation time histories in response 
to increasing magnitudes of crosswind gust ( 5 , 
10, 15, 20, 25, and 30 ft/s). 

So the shape-change effector array suite 
examined in this report is of limited practicality 
in the regard that it provides relatively low 
disturbance-rejection capability. Although the 
effector suite appears suitable to stabilize the 
lateral-directional dynamics in the presence of 
mild turbulence, and to perform relatively low- 
rate rolling maneuvers, the effector suite lacks 
sufficient authority to act as the sole means of 
controlling the vehicle in all situations. Rather, 
such a complement of shape-change effector 
arrays may be use to augment conventional 
devices, potentially providing a low-observable 
flight mode that may be flown through an 
autopilot executing only low-rate maneuvers. 
They may also be used to provide a certain 
degree of control redundancy, reconfigurability 
or damage tolerance as an augmentation to a 
conventional effector complement. There is also 
potential for new novel effector concepts to be 
developed that offer much greater authority 


through separation control, forebody vortex 
manipulation, or other means, and in such cases 
additional control schemes for employing these 
devices in the portions of the flight regime at 
which they are most useful will need to be 
devised. 


CQ NCLVPING REMARK S 

This investigation has examined the use of 
distributed shape-change effector arrays in a 
flight control system for an example aircraft 
design. The subject aircraft was Lockheed- 
Martin's Innovative Control Effector (ICE) 
configuration. The distributed shape-change 
effector arrays were modeled as a series of 
bumps normal to the aircraft surface which could 
be deployed to generate control moments. 

A sensitivity analysis was performed which 
provided insight regarding favorable locations 
for the placement of the shape-change device 
arrays on the ICE geometry. These sensitivity 
data were incorporated into an interactive shape- 
change array design tool, and a suite of arrays 
was developed and used in a lateral-directional 
control law design to illustrate allocation and 
deployment methodologies. A “quantized 
spanwise” deployment scheme was devised for 
activation of the device arrays. That is to say that 
each device in an array was either completely on 
or off (rather than proportionally deployed), and 
more of them were turned on to produce larger 
forces as needed, from the most inboard to the 
most outboard span station of each array. 

The control design was used in a dynamic 
simulation of the ICE aircraft for a flight 
condition of Mach 0.6 at 15,000 ft. Though the 
authority of the particular device arrays that were 
designed in this investigation was quite low, the 
simulation results indicate that the effector suite 
possessed sufficient authority to stabilize and 
maneuver the example vehicle model, executing 
relatively low-rate rolling maneuvers at 10 deg/s. 
When substantial atmospheric disturbances were 
included in the simulation (turbulence levels of 6 
ft/sec or cross wind gusts of >28 ft/sec), the 
device arrays did not possess sufficient authority 
to maintain stability of the vehicle's lateral- 
directional dynamics using the control law 
presented in this report. 

The shape-change device arrays studied in 
this investigation present a promising alternative 
for use in aircraft flight control. However, it is 
clear that experimental demonstrations of their 
application and further development of their 
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control authority would bolster their case. 
Future research will focus on experimental 
validation of the predicted authority of various 
flow control devices and on flight control using 
large arrays of interacting effector devices. A 
broader database that will characterize the 
variation of control authority of the shape- 
change device arrays with angle of attack and 
Mach number is also planned. 
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